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We study the relation between short-time vibrational modes and long-time relaxational dynamics 
in a kinetically constrained lattice gas with harmonic interactions between neighbouring particles. 
We find a correlation between the location of the low (high) frequency vibrational modes and regions 
of high (low) propensity for motion. This is similar to what was observed in continuous force 
systems, but our interpretation is different: in our case relaxation is due to localised excitations 
which propagate through the system; these localised excitations act as background disorder for the 
elastic network, giving rise to anomalous vibrational modes. Our results show that a correlation 
between spatially extended low frequency modes and high propensity regions does not imply that 
relaxational dynamics originates in extended soft modes. We consider other measures of elastic 
heterogeneity, such as non-affine displacement fields and mode localisation lengths, and discuss 
implications of our results to interpretations of dynamic heterogeneity more generally. 



Two of the central features of glass forming systems 
are, on the one hand, the increasingly fluctuating, het- 
erogenous dynamics [1 that accompanies the growth of 
relaxational timescales in supercooled liquids [2 , and on 
the other, the anomalies in their vibrational spectrum 
when considered as amorphous solids [3J H]. While the 
former pertains to structural relaxation on long times 
(say of the order of seconds for a liquid near the exper- 
imental glass transition temperature), and the latter to 
vibrational motion at high frequency (say in the THz for 
the Boson peak [3] of glasses) , an important question is 
whether these two characteristic aspects of the dynamics 
of glass formers are related, and in particular whether 
there is a common structural origin for both of them. 

Interesting recent studies of systems of hard spheres 
[5], hard disks [Bj, and soft disks |7J have revealed that 
the spatial localization of the anomalous low frequency 
modes of vibrations around either coarse-grained con- 
figurations in meta-basins [5J [B], or around inherent- 
structures [7J , give a good indication of the spatial distri- 
bution of structural relaxation at much longer times. It 
has been argued [5] [BJ that this correspondence is causal 
[7J, in that it is these soft vibrational modes which pro- 
vide the underlying structural mechanism to long time 
relaxation. 

In this paper we address this problem by studying 
spatial correlation between structural relaxation and vi- 
brational modes in suitably generalized kinetically con- 
strained models of glasses [5]. Our main results are il- 
lustrated in Fig. 1. The top left panel is a "propensity" 
map [9] [13] of a two- vacancy assisted lattice gas, or (2)- 
TLG [10], showing the average spatial distribution of 
structural relaxation starting from a typical configura- 
tion (details are given below). The right hand panels 
show the location of the lowest frequency modes (top) 
and highest frequency modes (bottom) of the elastic net- 
work obtained from the same configuration when neigh- 
bouring particles interact via harmonic springs. There is 
a clear spatial correlation between low frequency modes 



and high propensity for relaxation, and high frequency 
modes and low propensity for relaxation, as in atom- 
istic models [5j [6j [7] . In our model, however, relaxation 
is via the propagation of localized clusters of vacancies 
[TT1 [T2l [TB] , and does not emanate from the soft vibra- 
tional modes. In fact, it is precisely the presence of these 
localized defects in the elastic network that gives rise to 
the anomalous soft modes. We conclude that soft vibra- 
tional modes do indeed provide a good indicator of long 




FIG. 1: Clockwise from top left: (a) Propensity map of the 
(2)-TLG at density p = 0.75 averaged over 100 trajectories. 
Black circles indicate the isoconfigurational average distance 
travelled after time t = 5 x 10 3 . (b) The average participa- 
tion, ViiiS) of the lowest frequency modes, u < 0.1, and (c) 
the highest frequency modes, ui > 2.4, of the central force net- 
work, (d) The connected correlation, C(t,u>), between isocon- 
figurational persistence, Pi°(t), and the participation, Vi{(jj), 
as defined in the text. From bottom to top the frequencies 
are lo = 0.02,0.1,1,2,2.3,2.4. 
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time dynamic heterogeneity, but this correlation does not 
imply they are the underlying mechanism for long time 
relaxational dynamics in glass formers in general. 

The models we study are a generalization of con- 
strained lattice gases, first introduced as models of glasses 
in Ref. [TJ]. In these models hard-core particles oc- 
cupy the vertices of a lattice, with single occupancy per 
site, and with no static interactions between them. For 
the dynamics, particles try to hop to neighboring sites, 
but the local hopping rates depend on the occupancy 
of surrounding sites, so as to mimic steric interactions. 
These models have the trivial thermodynamics of a non- 
interacting lattice gas, but their dynamics at high densi- 
ties displays many of the features of glass formers, such 
as non-exponential relaxation jTTJ [T3], dynamic hetero- 
geneity , transport decoupling jH] , and aging [T5] (see 
Ref. [5j for a review) . 

In particular, we focus on the two-vacancy assisted lat- 
tice gas model on the triangular lattice [10], or (2)-TLG, 
where the kinetic constraint is explicitly due to steric re- 
strictions: a particle can hop to an empty nearest neigh- 
bour site only if the common two neighbouring sites are 
also empty. We also consider its three dimensional vari- 
ant on an FCC lattice, where the constraint is that the 
four common neighbours of the sites undergoing the tran- 
sition arc empty. We call this model (4)-FLG [16 . 

Figure 2 shows propensity maps for the (2)-TLG, 
which illustrate the spatial localization of relaxational 
dynamics in the model. The leftmost panel is a typical 
configuration of the model at a density, p — 0.75 in this 
case, for which relaxation is highly non-exponential and 
heterogeneous [IS]. Panels two to five show the average 
particle displacement at increasing times in the so-called 
iso-configurational ensemble [9 , i.e., for all trajectories 
starting from the same initial configuration (the one in 
the leftmost panel). For a detailed study of propensity 
in the (2)-TLG model see Ref. [13] . 

In order to study vibrations together with dynam- 




ical heterogeneity we generalize the models above by 
adding harmonic interactions between neighboring par- 
ticles. That is, any two occupied nearest neighbor sites 
interact through a linear spring of force constant k and of 
rest length equal to the lattice spacing. The vibrational 
Hamiltonian then reads, 



x 1 



rijk [(Sri - 8fj) ■ f° 



(1) 



where (ij) means that the sum is over nearest neighbor 
pairs, rij = 0, 1 indicates whether lattice site i is empty 
or occupied, <5r*i is the displacement of the particle whose 
equilibrium position is site i, and r?- is the unit vector be- 
tween sites i and j. Each configuration thus gives rise to 
a disordered elastic network due to the presence of vacan- 
cies in the particle configuration. This elastic problem is 
well known, as it corresponds to the "central force" prob- 
lem of rigidity percolation [17]. In particular, Eq. for 
the (2)-TLG and (4)-FLG models correspond to the site 
diluted triangular and FCC central force networks stud- 
ied in Ref. [18]. 

Figure 3 shows the average density of states (DoS), 
D(oj), of vibrational modes for elastic networks corre- 
sponding to the (2)-TLG model and (4)-FLG model at 
various densities. For densities lower than that of the 
full lattice, p < 1, the DoS presents an excess of low fre- 
quency modes. Debye scaling, D(u) oc a;" -1 , is recovered 
at low enough frequencies. The presence of excess modes 
becomes more pronounced as the density decreases. At 
a critical density p c the excess tail extends all the way to 
u> = 0. This is the isostatic point at which the system is 
marginally rigid 1 17] . For densities below p c there is an 
extensive number of zero frequency modes. We find that 
for the (2)-TLG p c « 0.7 and for the (4)-FLG p c » 0.5, 
which agrees with the results of Ref. [TH] for the triangu- 
lar and FCC lattices, respectively. 
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FIG. 2: Propensity maps for a (2)-TLG at density p = 0.75 
(second to fifth panels) for the same configuration as Fig. [I] 
(shown on the leftmost panel) averaged over 100 trajectories. 
Black circles indicate the average distance travelled by each 
particle after time, from left to right, t — 5 x 10 1 , 5 x 10 2 , 5 x 
10 3 , 5 x 10 4 . At this density, r Q « 10 4 . 



FIG. 3: Density of states for the central force networks of 
the (2)-TLG model (left panel) and the (4)-FLG model (right 
panel), at various densities p above the isostatic density p c . 
The curves at p = 1 are exact. The curves for all other 
densities were obtained by numerical diagonalization of the 
dynamic matrix of systems with 4000 normal modes. 
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For all densities below that of the full lattice, p < 1, 
the vibrational modes are spatially distributed in a non- 
trivial manner. Figure 1 illustrates this. It shows the 
spatial weight, Vi(uj), of modes of low frequency (Fig. 1, 
top-right) and high frequency (Fig. 1, bottom-right) for 
the (2)-TLG configuration of Fig. 2, withfj(w) = e^-e^ 
where is the two dimensional i-th component of 
the eigenvector of frequency to. With this definition 
J2i — 1 for all to. Fig. 1 (top-right) also shows 

the propensity map for the same configuration (cf. Fig. 
2) overlayed on the low frequency eigenvectors. The pic- 
ture suggests a close correlation between areas of high 
propensity for relaxational motion and the location of 
soft modes. Fig. 1 (bottom-right) shows a similar spa- 
tial correlation between regions of low propensity (shown 
as the inverse of Fig. 1, top- left) and the location high 
frequency modes. 

In order to quantify this spatial correlation we define 
the following cross correlation between local propensity 
for motion and vibrations: 

C(t,u;) = (pl c (t)^))-(p] c (t)). (2) 

Here pi (t) denotes the persistence function [12] of particle 
i, i.e. Pi(t) — 1 if particle i has not moved up to time t, 
and Pi(t) = otherwise. p\ (t) is the iso-configurational 
average of Pi(t), i.e. the average of the persistence field 
over all trajectories that start from a given configuration 
1 1 3 j - The average (•) in Eq. ((|) is over all equilibrium 
configurations at a given density. Fig. 1 (bottom-left) 
shows C{t,u}) as a function of time t for various vibra- 
tional frequencies for the (2)-TLG at density p — 0.75. 
For high frequencies the correlation is positive, indicat- 
ing that particles which are more persistent than aver- 
age also participate in high frequency vibrational modes. 
For low frequencies the correlation is negative, indicating 
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FIG. 4: Peak in the absolute value of the correlation C(t,ui) 
between propensity for relaxational motion and vibrational 
modes, in the (2)-TLG and (4)-FLG, as a function of den- 
sity; cf. Fig. 1 (bottom-left). The correlation between fast 
regions/high-frequency modes decreases with increasing den- 
sity. The correlation between slow regions/low-frequency 
modes increases with increasing density. 



that fast relaxing particles (low p\ c ) are also those which 
participate in soft vibrations. In both cases the correla- 
tion is non-monotonic in time, peaking at times around 

Figure 4 shows how the spatial correlation between vi- 
brations and dynamic heterogeneity depends on density. 
It plots the peak value of C(t 7 uj) for vibrational modes 
of the lowest and highest non-trivial frequencies acces- 
sible in the simulations. The fast-relaxation/soft-mode 
correlation increases with density for all densities larger 
than the isostatic one, p > p c . In contrast, the slow- 
relaxation/high- frequency correlation decreases with in- 
creasing density. These trends are similar in dimension 
two, (2)-TLG, and dimension three, (4)-FLG. Below the 
isostatic point the correlation between slow regions and 
low frequency modes changes sign as the vibrational spec- 
trum becomes plagued by zero modes. 

The heterogeneity of vibrational modes in the (2)-TLG 
and (4)-FLG increases with decreasing particle density. 
This vibrational heterogeneity can be seen by computing 
the localization length of each mode. We do so by us- 
ing the method of Ref. [19 , where a mode localization 
length is extracted by tracking the change in its eigen- 
value due to asymmetric perturbations of the dynamical 
matrix |20j . Figure 5 (top) shows the change in the local- 
ization length computed in this way for different densities 
in the (2)-TLG. In the full lattice limit, p—1, there are 
no vibrational anomalies and all the modes are extended. 
As soon as a small density of vacancies is present some 
of the modes become localized, as illustrated in Fig. 5 for 
p = 0.95. More modes become localized as the density 
of particles decreases. Near to the isostatic point most 
modes have localization lengths that are smaller than our 
system size, see middle panel of Fig. 5 (top) . Beyond the 
isostatic density the system becomes elastically unstable. 

Elastic heterogeneity also becomes evident in the non- 
affine response of our system to an external strain 
[2T1 122] . Following Ref. [25], we deform the system 
through a uniform strain in the x-direction, e <C 1, 
by rescaling all particle positions in an affine manner, 
r~io —* r* i0 . After the stretch we minimize the energy us- 
ing conjugate gradient descent. Because of the periodic 
boundary conditions, once the affine displacement has 
been made the system will not be able to relax back to 
its original configuration. The final positions of the par- 
ticles can be expressed as the sum of an affine part, r[ 
and a nonaffine part, u,. 

For high densities the conjugate gradient algorithm is 
very efficient at finding the minimum energy configura- 
tion. As density is lowered towards the rigidity threshold 
the time to converge starts to increase rapidly. This is 
due to the formation of weakly connected clusters that 
can make large displacements at a small energy cost. 
These clusters are the beginnings of the floppy modes 
that are able to operate almost independently from the 
rest of the system. At the isostatic point gradient de- 
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scent cannot converge. To prevent these divergences we 
have added a small confining harmonic potential of force 
constant at each site to Eq. ([IJ . This extra term is di- 
agonal, and so the normal modes remain unchanged but 
with a shift in frequency. 

Figure 5 (bottom) shows the average magnitude of the 
non-affine response, ((ui/e) 2 ) 1 / 2 , as afunction of the den- 
sity in the (2)-TLG model. Close to the isostatic point 
the average non-affine deformation increases very rapidly, 
and would appear to diverge when f2 — > 0. The inset to 
Fig. 5 (bottom) has the projection of the non-affine de- 
formation on the normal modes, (u\u>) = X^Uj ■ e, w , as 
a function of frequency of the modes, for the (2)-TLG 
at density p — 0.75, showing that non-affinity is carried 
preferably by softer modes. 

In summary, we have studied the correlation between 
dynamic heterogeneity and anomalous vibrations in a two 
and a three dimensional constrained lattice gas model 
of glasses. The structural relaxation of these models at 
high density is similar to that of glass formers, display- 
ing non-exponential relaxation and dynamic heterogene- 
ity [HI [13l [16] . Their vibrational properties are those of 
well-studied random networks |18j . and mimic character- 
istic aspects of the anomalous vibrations of glasses: ex- 
cess low frequency modes, non-affinity and elastic hetero- 
geneity, all related to the presence of an isostatic point. 
We have found that the location of anomalous vibrational 
modes correlates to dynamic heterogeneity of structural 
relaxation, as is observed in atomistic systems [5JIH1[7]. 
In our case, however, structural relaxation, and there- 
fore dynamic heterogeneity, originate in localized vacan- 
cies [HJII3], and not in the extended structures that the 
soft-modes span. In fact, vacancies act as quenched lo- 
calized defects for the vibrations, cf. Eq. (JT|), giving rise 
to the anomalous elastic behaviour observed. We have 
thus shown through these simple examples that a corre- 
lation between soft modes and propensity does not imply 
a causal relation for relaxation mechanisms. A similar 
situation may hold in atomistic models as well. 

This work was supported by EPSRC grant 
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